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^ Abstract 



We present a nonperturbative dynamical study of production in proton- 



nucleus and nucleus-nucleus collisions at SPS energies on the basis of the covariant 
transport approach HSD. For p + W reactions at 200 GeV bombarding energy the 
dimuon yield for invariant masses m < 1.6 GeV is found to be dominated by the 
decays of the T],p,u and mesons. For 200 GeV/A S + W collisions, however, the 
dimuon yield shows an additional large contribution from tt + tt~, K + K~ and 7rp 
channels. We find that for 'free' meson masses and form factors the experimental 
cross section is clearly underestimated for S + W in the invariant mass range 0.35 
GeV < m < 0.65 GeV and that the HELIOS-3 data can only be reproduced within 
a hadronic scenario, if the p/w-meson mass drops with baryon density. This finding 
suggests a partial restoration of chiral symmetry in S + W collisions at SPS energies. 
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The question of chiral symmetry restoration at high baryon density is of fundamental 
interest since a couple of years |l|, 0, but a clear experimental evidence has not been 
achieved, so far. The enhancement of strangeness production as e.g. seen in the AGS 
data for the K + /tt + ratio || might be a signiture for such a transition EL however, other 
hadronic scenarios can be cooked up to describe this phenomenon as well || . On the other 
hand, dileptons are particularly well suited for an investigation of the violent phases of 
a high-energy heavy-ion collision because they can leave the reaction volume essentially 
undistorted by final-state interactions. Indeed, dileptons from heavy-ion collisions have 
been observed by the DLS collaboration at the BEVALAC || ^, |§ and by the CERES 
IJ and HELIOS collaboration [10, 11 at SPS energies. 

The enhancement of the low mass dimuon yield in S + W compared to p + W col- 
lisions [ID] has been first suggested by Koch et al. [T2j to be due to 7t + tt~ annihilation. 
Furthermore, Li et al. E3[ have proposed that the enhancement of the e + e~ yield in S + 



Au collision as observed by the CERES collaboration || should be due to an enhanced 
p-meson production (via 7r + 7r~ annihilation) and a dropping p-mass in the medium. In 
fact, their analysis - which was based on an expanding fireball scenario in chemical equi- 
librium - could be confirmed within the microscopic transport calculations in ref. JTjJ] ; 
however, also a more conventional approach including the increase of the p-meson width 
in the medium due to the coupling of the p, it, A and nucleon dynamics ]15[ 16], [17| was 
found to be compatible with the CERES data. In this paper, we will carry out a similar 
study as in ref. [14]] but for the dimuon data of the HELIOS-3 collaboration [JTOJ which 
provide independent information on dilepton production with rather good statistics. 

In continuation of the dilepton production studies in refs. |14|, [18], [19], ^(J we describe 
the dynamics of proton-nucleus or nucleus-nucleus reactions by a coupled set of covariant 
transport equations with scalar and vector self-energies of all hadrons involved. Explicitly 
propagated are nucleons, A's, N*(1440), N*(1535) resonances as well as 7r's, r/s, p's, u;'s, 
$'s, kaons and K*'s with their isospin degrees of freedom. For more detailed information 
on the self-energies employed we refer the reader to ref. where the transport approach 
HSD[] is formulated and applied to nucleus-nucleus collisions from SIS to SPS energies. 

In this analysis we calculate p + p~ production taking into account the contributions 
from nucleon-nucleon, pion-nucleon and pion-pion bremsstrahlung, the Dalitz-decays r] — > 
7P + p~ and uj — > 7r°/i + p~, the direct dimuon decays of the vector mesons p, u, $ as well 
as 7r + 7r~, K + K~ and 7rp production channels. The nucleon-nucleon, pion-nucleon and tttt 

1 Hadron-String-Dynamics 



bremsstrahlung are evaluated in a phase-space corrected soft photon approximation as 
described in refs. [18, SO]. Though this limit is questionable for the bombarding energy 



of interest, the soft photon approximation provides an upper limit for the bremsstrahlung 
contribution |^TJ, which will turn out to be small in comparison to the meson decay 



channels. 

The pion annihilation proceeds through the p-meson which decays into a virtual mas- 



sive photon by vector meson dominance. The cross section is parametrized as [18, pQ, 22 



S \mJ V m 

where the free form factor of the pion is approximated by 

\FAm)\ 2 = — 2 / . (2) 

[m A — m p 1 ) 1 + rripl z 

In eq. (2) m is the dilepton invariant mass, a is the fine structure constant, and 

m p = 775 MeV , T p = HSMeV . 

So far, eq.(|I|) describes the free pion annihilation cross section; possible medium modifi- 
cations of the p-meson will be discussed below. 



The cross section for the K + K annihilation is parametrized as p3 



a ^-^- (m) = ^(-) 2 J 1 -^ lFK{m)P , (3) 
S \mj V mr 

where the free form factor of the kaon is approximated by 

\Mm)\ a = ~ TL 2r2 ( 4 ) 

9 (m 2 - m%y + m%Y% 

with 

m$ = 1020Me1/ , T$ = AA3MeV . 

The cross section for dimuon production in vr + p~, 7r~p + scattering is given by the 
Breit-Wigner cross section [2~3| for the formation of a $-meson times the partial decay of 
the $ into a p + p— pair, 
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with = 0.13 and -B OU £ = 2.5 xl0~ 4 . We note that the resonance approximation (5) 
also works very well for K + K~ annihilation as described by eq. (3) when including a 
further spin factor of 3 and B in = 0.49. 

The uj Dalitz-decay, furthermore, is given by |25| : 




where the form factor is parametrized as 



A1 



m 
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with 



A, = 0.60 GeV , I\ = 75 MeV. 



The form (7) for F u is actually fitted to the dimuon data for p + W at 200 GeV for 
0.4 GeV < m < 0.6 GeV and slightly deviates from the parametrization in |^3] (with 
A s = 0.65 ± 0.03 GeV) but is still compatible with the data presented in [25|. We note 
that this formfactor is of interest in its own since it shows a significant enhancement 



compared to the vector-meson-dominance model p5[ . More accurate experimental data 
on this formfactor as well as further theoretical studies will be needed to clarify the 
underlying physics. 

The direct decays of the vector mesons to fi + fi~ are taken as 
da , , 1 myTy T v _ 
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with r*_ M+M -/r^ = 2.5xio- 4 , r^ +p -/r^ = 7.1 xio- 5 and iy^ +M -/rj 

x 10 -5 . Though of minor importance we also include the direct decay of the r/- meson with 
r ?? _ >/Lt + jU - /r*°* = 5.7 xlO -6 . The charged particle multiplicity N c - which also provides 
the normalization of the experimental data - is recorded in the pseudorapidity interval 
3.7 < r] < 5.2. 

Transport models similar to ours have been applied for e + e~ production at BEVALAC 

I] and at AGS energies in [R0]. The first full microscopic 



energies in refs. [26, 27, 28 



studies for e + e~ production at SPS energies have been reported in [14|; the same covariant 
transport approach is also used in our present analysis. 

We have calculated the yield for p + W at 200 GeV using the experimental cuts 

in rapidity y and transverse mass rriT of the dimuon source, 



m T >4(7-2y), rrvr > A (2m„)» + ( ^^ £ ) \ ( 10 ) 

as well as the experimental resolution in energy [TO . 

Our results are displayed in Fig. 1 for p + W at 200 GeV in comparison with the 
HELIOS-3 data fpL0| . In the upper part of Fig. 1 we show the calculated results with 
an energy resolution of Am = 15 MeV, while the experimental resolution was used in 
the lower part of Fig. 1 in order to demonstrate the sensitivity of the spectra to the 
actual mass resolution in the experimental data. The full solid curves in Fig. 1 display 
the sum of all individual contributions which is dominated by the decays of the mesons. 
The bremsstrahlung contributions (7rN, pN and nn channels summed up in the single line 
denoted by 'brents') as well as 7r + 7r~ annihilation (short dashed line) are of minor impor- 
tance for the proton induced reaction in full analogy to the CERES e + e~ data |5|, TA\. The 
very good reproduction of the HELIOS-3 data, where medium effects are negligible, allows 
us to perform a related analysis of nucleus-nucleus collisions with particular emphasis on 
medium properties of the reaction channels involved. 

We thus go over to the system S + W at 200 GeV/A where we expect to describe 
the global reaction dynamics with the same quality as for S + Au at 200 GeV/A (cf. 
Fig. 2 of ll^]). The results of our calculation, where no medium effects are incorporated 
for all mesons, are displayed in Fig. 2 (thick solid line) in comparison to the data [|i~0|j . 
The individual contributions from bremsstrahlung channels (denoted by 'brems'), r\ and 
uj Dalitz decays as well as direct vector meson decays are explicitely indicated in Fig. 2. 
Concentrating on the invariant mass range m < 1 GeV we find the experimental yield to be 
underestimated for 0.35 GeV < m < 0.65 GeV and slightly overestimated for m pa 0.8 GeV. 
In the $-mass region the cross section includes the direct $ production channels as well 
as those from and K + K~ collisions^]. The tt + tc~ annihilation component (short dashed 
line) is found to be of the same order of magnitude as the direct p-decay (long dashed 
line) within the experimental cuts and significantly smaller than in case of the CERES 
data |14| - which focus on the mid-rapidity regime - while the additional bremsstrahlung 
contributions again are of minor importance. 



2 The kaon production channels irir — > KK are included in the transport calculation 



Since the reproduction of the HELIOS-3 data in Fig. 2 is rather poor we will now 



explore possible in-medium modifications of the mesons as advocated in refs. [14. |2(]. 23 
|3"T| ]. Before going over to the actual model, we show in Fig. 3 for S + W at 200 
GeV/A the p-mass distribution versus the baryon density p# in units of po ~ 0.168/m~ 3 , 
where is baryon density is computed in the local rest frame, i.e. Pb{x) = jfxj^, while 
jn{x) denotes the Lorentz invariant baryon current. In Fig. 3 the solid line represents 
the distribution of p's produced in baryon-baryon and pion-baryon collisions whereas the 
dashed line corresponds to those from 7r + 7r~ annihilation. We note that the distributions 
in Fig. 3 are recorded at the creation time of the p-mesons, whereas the distribution 
at the actual decay time is slightly shifted to smaller density. Furthermore, only those 
p-mesons have been considered which fulfill the experimental cuts (10) in order to provide 
a close connection to the HELIOS-3 experiment. These distributions clearly indicate that 
most of the p's are produced at nonzero baryon density and even for ps > 3po- 

We now turn to the possible modifications of the p-meson in the medium. From QCD 
inspired models [[[], |31| or estimates based on QCD sum rules |32], |33|] it has been predicted 
that the p-meson mass decreases with density. In order to explore the compatibility of 
such scenarios with the HELIOS-3 data, we have performed calculations with a medium- 
dependent p- (and c<j)-mass according to Hatsuda and Lee [ 3"2j , i.e. 

m* « m° p (1 - 0.18p B /p ) > m u + m d « UMeV, (11) 

where ps(t) is the actual baryon density during the decay of the p-meson. A similar 



scaling with density has also been predicted by Shakin et al. |3"I| . Since the width of 
the p-meson is determined by the decay into two pions, which are propagated without 
selfenergies as in [f|], we have parametrized the in-medium width T* as 

YYl*2 — ATP? 

r ^ r '7^4 e ^~ 2m ^ (12 > 

The dropping of the p-meson mass is associated with a scalar self-energy of the meson 
which is determined by the local baryon density; the propagation of the 'quasi-particle' 
with effective mass m* thus couples to the baryon current during the expansion, and the 
meson becomes 'on-shell' asymptotically due to a feedback of energy from the mean fields, 
which thus ensures that the total energy of the system is conserved. 

The results of this simulation are shown in Fig. 4 in comparison with the HELIOS-3 
data. The ir + 7T~ annihilation component (short dashed line) now shows a peak at invariant 
masses m k 0.6 GeV and leads to a significant enhancement of low mass dileptons. A 



further enhancement at lower masses is also provided by the direct p-decays (long dashed 
line) due to their shifted masses. In comparison to Fig. 2 the experimental dimuon 
spectrum is now much better reproduced for invariant masses m < 1 GeV; the missing 
yield for m > 1.2 GeV might be due to direct charm production channels or higher meson 
resonances which are not explicitly included in the present calculations. 

Of further interest is the experimental rapidity dependence of the dimuon spectra. 
Here, explicit data are available within the cuts (10) for y < 3.9, 3.9 < y < 4.4, and y > 
4.4 . The comparison of our calculations with the data is shown in Fig. 5 for the rapidity 
cuts 3.9 < y < 4.4 and y > 4.4. The dotted lines correspond to calculations with free 
meson masses and formfactors whereas the solid lines are the result within the dropping 
meson mass scenario (11) with the in-medium width (12). Again the computations with 
free meson masses fail to reproduce the data especially for the lower rapidity interval 
from 0.35 - 0.65 GeV and about 0.8 GeV whereas the agreement is remarkably fine when 
including the shift of the meson masses with density. 



As advocated by Haglin [34| the width of the p-meson should again increase in a 
dense mesonic environment due to the mesonic interaction channels. In order to explore 
such effects we have performed calculations with a dropping p-mass according to (11), 
however, keeping its width at the free value, i.e. T* = T p . The results of this 'theoretical 
experiment' are displayed in Fig. 5 in terms of the dashed lines. Whereas for y > 4.4 the 
agreement with the data improves, the yield is slightly underestimated for m ^ 0.5 

GeV in the lower rapidity interval. We thus cannot extract a clear evidence for a reduced 
width of the p-meson at high baryon density from the present experimental data. 

Within the same spirit we have also investigated, if the broadening of the p-meson - 



without a shift of its pole - according to Herrmann et al. ||15| , which is due to the coupled 
7r and A-hole dynamics, might improve the description of the data. However, contrary 
to the e + e~ spectra of the CERES-collaboration (cf. ref. [[H]), there is no significant 
enhancement of the spectra for 0.3 GeV < m < 0.6 GeV (in comparison to Fig. 2) 
such that this mechanism can no longer be considered as dominant. Furthermore, since 
the pion densities in the rapidity regime of interest are large compared to the baryon 
densities (cf. Fig. 2 of |nj), the pion selfenergy effects from A-hole loops are most likely 
overestimated |B5|. 

In summary, we have studied p + p~ production in proton and heavy-ion induced re- 
actions at 200 GeV/A on the basis of the covariant transport approach HSD Q. We 
have incorporated the contributions from proton-nucleon, pion-nucleon and pion-pion 
bremsstrahlung, the Dalitz-decay of the A, 77 and u as well as 7r + 7r~ annihilation and the 



direct dilepton decay of the vector mesons p,u, $ as well as K + K~ and np channels. It 
is found that for p + W at 200 GeV the mesonic decays almost completely determine 
the dilepton yield, whereas in S + W reactions the it + tt~, K + K~ annihilation channels 
and irp collisions contribute substantially The experimental data taken by the HELIOS-3 
collaboration jnj are underestimated by the calculations for invariant masses 0.35 GeV 
< m < 0.65 GeV when using free form factors for the pion and p- meson. 

We have, furthermore, examined if a shift of the p- (and u;)-mass according to QCD 
sum rules - as suggested by Hatsuda and Lee [0 or Shakin et al. |31] - is compatible 



with the data as in case of the e + e~ spectra of the CERES collaboration jyj]. Indeed the 
description of the p + p~ spectrum improves substantially within this 'dropping mass' sce- 
nario at all rapidity bins. However, the effects at forward rapidities are not as pronounced 
as close to midrapidity where the dilepton yield is clearly dominated by the 7r + 7r~ anni- 
hilation channel (cf. Figs. 3-5 of [0). Nevertheless, the explicit rapidity dependence of 
the experimental dimuon spectra, which is in line with the dropping meson mass scenario, 
provides a strong hint for the partial restoration of chiral symmetry at high baryon den- 
sity. Data with high statistics especially at midrapidity or alternatively central collisions 
of Pb + Pb at SPS-energies are expected to allow for more definite conclusions about the 
chiral phase transition of hadronic matter at 4-6 xpo baryon density. 
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Figure Captions 

Fig. 1: Comparison of our calculations for the differential dimuon spectra (thick solid 
line) with the experimental data [I(J for p + W at 200 GeV. Upper part: with an energy 
resolution of 15 MeV; lower part: with the experimental resolution. The individual con- 
tributions from the vector mesons and Dalitz decays as well as bremsstrahlung channels 
are directly indicated in the figure. The 7T + ir~ annihilation channel and the decay from 
p-mesons produced in baryon-baryon and meson-baryon collisions are shown by the short 
dashed and long dashed lines, respectively. 



Fig. 2: Comparison of our calculations for the differential dimuon spectra (upper thick 
solid line) with the experimental data JTO for S + W at 200 GeV/A when employing no 



medium modification of the mesons. The 7r + 7r~ annihilation channel and the decay from 
p-mesons produced in baryon-baryon and meson-baryon collisions are shown by the short 
dashed and long dashed lines, respectively. 



Fig. 3: The p-mass distribution as a function of the baryon density p# (in units of po) for 
a central S + W reaction at 200 GeV/A within the HELIOS-3 acceptance (10). Solid line: 
p-mesons from baryon-baryon and meson-baryon reaction channels; dashed line: p-mesons 
from 7r + 7r~ annihilation. 



Fig. 4: Comparison of our calculations for the differential dimuon spectra (thick solid 
line) with the experimental data [[K]] for S + W at 200 GeV/A when including a shift 
of the p-meson mass (11) according to the prediction by Hatsuda and Lee and a 
decrease of its width according to (12). The 7r + 7r~ annihilation channel and the decay 
from p-mesons produced in baryon-baryon and meson-baryon collisions are shown by the 
short dashed and long dashed lines, respectively. 



Fig. 5: Result of our calculations for the differential dimuon spectra for the rapidity cuts 
3.9 < y < 4.4 and y > 4.4 for S + W at 200 GeV/A in comparison to the HELIOS-3 data 



1C]. The dotted lines correspond to a calculation with free meson masses and formfactors, 



the solid lines to a calculation with a dropping p-mass according to (11) and width (12), 
whereas the dashed line is obtained in the dropping p-mass scenario (11) with the free 
p-meson width. 
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